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The spin-1/2 Ising-Heisenberg model with the pair XYZ Heisenberg interaction and quartic Ising 
interactions is exactly solved by establishing a precise mapping relationship with the corresponding 
zero-field (symmetric) eight-vertex model. It is shown that the Ising-Heisenberg model with the 
ferromagnetic Heisenberg interaction exhibits a striking critical behavior, which manifests itself 
through re-entrant phase transitions as well as continuously varying critical exponents. The changes 
of critical exponents are in accordance with the weak universality hypothesis in spite of a peculiar 
singular behavior to emerge at a quantum critical point of the infinite order, which occurs at the 
isotropic limit of the Heisenberg interaction. On the other hand, the Ising-Heisenberg model with 
the antiferromagnetic Heisenberg interaction surprisingly exhibits less significant changes of both 
critical temperatures as well as critical exponents upon varying a strength of the exchange anisotropy 
in the Heisenberg interaction. 
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I. INTRODUCTION 

Over the last six decades, exactly solvable models 
of interacting many-body systems have attracted con- 
siderable attention due to their uncoverable role in a 
fundamental understanding of order-disorder phenomena 
[ill 0i Q- particular, exactly solvable quantum spin 
models are currently at a forefront of theoretical research 
interest as they often exhibit a peculiar interplay between 
quantum fluctuations and cooperativity 0) S 0] ■ Despite 
extensive studies, however, phase transitions and criti- 
cal phenomena of rigorously solved quantum spin mod- 
els still belong to the most challenging unresolved issues 
to deal with, since quantum fluctuations usually prefer 
a lack of spontaneous order. Indeed, there are only few 
exactly solved quantum spin models that simultaneously 
exhibit both spontaneous long-range order, as well as, 
obvious quantum manifestations. 

It has been recently demonstrated that the sponta- 
neous order, which is accompanied with obvious quan- 
tum manifestations, might be an inherent feature of the 
hybrid Ising-Heisenberg planar models [13] whose lattice 



Note that all aforementioned Ising-Heisenberg planar 
models have exactly been solved with the help of suit- 
able algebraic transformations such as the decoration- 
iteration [l3] or the star-triangle |18| transformation, 
which establish a precise mapping equivalence with the 
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corresponding spin-1/2 Ising model. It is worthwhile to 
recall, moreover, that the decoration-iteration and star- 
triangle mapping transformations were substantially gen- 
eralized by Fisher [l9lj (see for details also Ref. ^2^) who 
firstly pointed out that in principle arbitrary quantum- 
mechanical system coupled to two or three outer Ising 
spins can be replaced via appropriate algebraic trans- 
formation by the equivalent expression containing pair- 
wise spin-spin interactions between the outer Ising spins. 
In such a way, one effectively establishes a precise map- 
ping relationship that connects the exact solution of some 
original model (which might describe a rather complex 
quantum-mechanical system) with the exact solution of 
the corresponding spin-1/2 Ising model, which is gener- 
ally known for many planar lattices of different topologies 
ri,"3, 20, 21, 22, 23]. 

On the contrary, there does not exist a general alge- 
braic transformation for any quantum-mechanical system 
coupled to four or more outer Ising spins if one consid- 
ers pairwise interactions between the outer Ising spins 
only [1^. Rojas et al. [2^ has recently found an exact 
evidence that it is nevertheless possible to include mul- 
tispin interactions between the outer Ising spins into the 
algebraic transformation in order to ensure its general 
validity. Thus, the algebraic transformation with effec- 
tive pair and quartic interactions is for instance required 
when searching for an exact treatment of the quantum- 
mechanical system coupled to four outer Ising spins in an 
absence of the external field. The spin-1/2 Ising model 
with pair and quartic interactions is however nothing 
but the alternative definition of the general eight-vertex 
model [25l . [2^ , which becomes exactly tractable by im- 
posing a special additional constraint to its vertex ener- 
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gies (Boltzmann's weights) known either as the Baxter's 




281 or the free- 
Bearing all this 



zero-field (symmetric) condition ll| 
fermion condition of Fan and Wu [2 
in mind, one could intuitively expect that there might ex- 
ist a certain class exactly solvable quantum-mechanical 
spin models for which the appropriate algebraic transfor- 
mation yields a precise mapping correspondence to the 
eight- vertex model generally satisfying the Baxter's zero- 
field condition P, [23, [2^ . The main goal of the present 
work is to show that the spin-1/2 Ising-Heisenberg model 
with the pair XYZ Heisenberg interaction and the quar- 
tic Ising interactions falls into this class of fully exactly 
soluble models. 



An importance of exact results to be obtained for the 
Ising-Heisenberg model with the pair and quartic inter- 
actions should be viewed in more respects. First, it is 
well known that the quartic interaction basically affects 
magnetic properties of several insulating cop per com- 
pounds such as La2Cu04 [3ll,[33,[3l,[33,|33,|3l| (undoped 
parent compound for high- Tc cuprates), SrCu203 [34| . 
La6Ca8Cu2404i [S^l, and La4SrioCu2404i [3l|. Even 
though it would be rather striking coincidence if some 
real magnetic material would obey very specific topo- 
logical requirements of the model under investigation, 
it is quite reasonable to suspect that our exact re- 
sults might at least shed light on some important as- 
pects of the critical behavior of real magnetic materi- 
als. Second, our model system represents a rare exam- 
ple of the exactly solved lattice-statistical model, which 
contradicts the standard universality hypothesis in that 
its critical exponents vary continuously with the inter- 
action parameters over the full range of possible val- 
ues of the critical exponents. Third, the exactly solu- 
ble vertex models have found over the last few decades 
manifold applications in seemingly diverse research ar- 
eas. Exact solutions of vertex-like models essentially 
tackles the problem of the residual entropy of two- 
dimensional ice [H, lioll ■ the Slater model of hydrogen- 
bonded ferroelectrics j4ll. [i^ \a§^ . the integrable quan- 
tum spin models such as the quantum Heisenberg chain 
m, m. El lipJilEiJii, the ice-type sohd-on-solid 
(SOS) modellHlla^^ the proble m of counting 



domino tilings 30, 3, 3, 3, , the three- and four- 
coloring problems of the square and hexagonal lattices 
[11, Us USUI, etc. 



The rest of this paper is organized as follows. In 
the next section, we will provide the detailed descrip- 
tion of the hybrid Ising-Heisenberg model and then, basic 
ideas of the exact mapping to the zero-field eight-vertex 
model will be explained. This is followed by the presenta- 
tion of the most interesting results for the ground-state 
and finite-temperature phase diagrams, which are sup- 
plemented by a detailed analysis of a rather strange non- 
universal behavior of critical exponents. Finally, some 
concluding remarks are given in Section [IVI 



II. ISING-HEISENBERG MODEL AND ITS 
EQUIVALENCE TO THE ZERO-FIELD 
EIGHT- VERTEX MODEL 

Consider a two-dimensional lattice of edge-sharing oc- 
tahedrons, each of them having four Ising spins cr = 1/2 
in a basal plane and two Heisenberg spins S = 1/2 in 
apical positions, as it is schematically depicted in Fig. [T] 
Suppose furthermore that each edge of the octahedron, 
which connects two Ising spins, is a common edge of two 
adjacent octahedrons (thin solid lines in Fig.[T]). The en- 
semble of all Ising spins then forms a square lattice, which 
has a couple of the Heisenberg spins above and below a 
center of each elementary square face formed by four Ising 
spins. Let both apical Heisenberg spins interact together 
via the pairwise XYZ Heisenberg interaction, while they 
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FIG. 1: The elementary unit cell of the spin-1/2 Ising- 
Heisenberg model. Full (empty) circles denote lattice posi- 
tions of the Heisenberg (Ising) spins, thick solid line repre- 
sents the pairwise XYZ Heisenberg interaction between the 
apical Heisenberg spins and both types of broken lines con- 
nect spins involved in the quartic Ising interactions. Thin 
solid lines connecting four Ising spins are guide for eyes only. 

are also assumed to be engaged in two different quartic 
Ising interactions with two Ising spins from opposite cor- 
ners of a square face (see Fig. [T]). The total Hamiltonian 
can be for convenience written as a sum over all elemen- 
tary unit cells (octahedrons) H = J2p'^p^ where each 

octahedron-cluster Hamiltonian Tip contains one pair in- 
teraction between the apical Heisenberg spins and two 
quartic interactions between the Heisenberg spins and 
their four Ising neighbors 
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(1) 



Above, the interaction parameters Jx, Jy, Jz denote spa- 
tial components of the anisotropic XYZ interaction be- 
tween the Heisenberg spins, while the interaction param- 
eters Ji and J2 label two quartic Ising interactions be- 
tween both apical Heisenberg spins and two Ising spins 
from opposite corners of a square face along two different 
diagonal directions (see Fig.[T|). 

It is of principal importance that the cluster Hamil- 
tonians of two different octahedrons commute with each 
other, i.e. [7ii,7ij] = 0, which immediately allows a par- 
tial factorization of the total partition function into a 
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product of cluster partition functions 

^ = EnTrpexp(-/37ip). (2) 

The summation X]{o-} emerge in Eq. Q is carried 
out over all possible configurations of the Ising spins, 
the symbol Tr^ denotes a trace over spin degrees of free- 
dom of the Heisenberg spin pair from the pth octahe- 
dron, P — l/(fcBT), /cb is Boltzmann's constant and T 
is the absolute temperature. After performing the rel- 



evant trace over spin degrees of freedom of the Heisen- 
berg spins, the partition function of the Ising- Heisenberg 
model can be rewritten into the form 

2^ = X^n^P ('^Pl'^p2)Crp3,Crp4) . (3) 
{a} P 

Apparently, the effective Boltzmann's factor ujp assigned 
to the pth octahedron now explicitly depends just on four 
Ising spins CTpi, CTp2, CTpS and (7p4 from its basal plane 
through the relation 



U!p{a, b, c, d) = 2 exp {Kz + Kiac + K2bd) cosh [Kx - Ky) + 2 exp [-Kz - Kiac - K2bd) cosh {K-j, + Ky) , (4) 



where we have introduced a new unified notation for the 
coupling constants Ka = PJa/'i (a = x, y, z, 1, 2) in or- 
der to write the Boltzmann's factor ([4]) in a more abbre- 
viated and elegant form. 

At this stage, the model under investigation can 
be rather straightforwardly mapped to the eight-vertex 
model on its dual square lattice. As a matter of fact, 
the product emerging in Eq. ([3]) can alternatively be 
performed over all elementary squares of the Ising spins 
forming a square lattice and the Boltzmann's factor (|4|) 
is invariant under the reversal of all four Ising spin vari- 
ables. In this respect, there are at the utmost eight differ- 
ent spin arrangements that have different energies (Boltz- 
mann's weights) and these can readily be related to the 
Boltzmann's weights of the eight-vertex model on a dual 
square lattice by the following procedure. If, and only if, 
the Ising spins located at adjacent corners of a square face 
are aligned opposite to each other, then solid lines are 
drawn on edges of a dual square lattice, otherwise they 
are drawn as broken lines. This is actually one of many 
alternative definitions of the eight-vertex model, since an 
even number of solid (broken) lines is always incident to 
each vertex of a dual square lattice. The diagrammatic 
representation of eight possible spin arrangements and 
their corresponding line coverings is shown in Fig. [21 It 
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FIG. 2: Eight possible Ising spin configurations and their 
relation to line coverings of the corresponding eight-vertex 
model on a dual square lattice. The sign '±' denotes the spin 
state cr^ = ±1/2. 

can be easily understood that each of eight possible line 
arrangements around a vertex of the dual lattice corre- 
sponds to two spin configurations, one is being obtained 
from the other by reversing all four Ising spins located at 



corners of an elementary square face. The hybrid Ising- 
Heisenberg model thus turns out to be equivalent with 
the eight-vertex model. With regard to this equivalence, 
the partition function of the Ising-Heisenberg model can 
be expressed in terms of the partition function of the 
eight-vertex model on a square lattice 

Z{T, Jx, Jy, Jz, Jl, Jl) — 2^8 — vertex 

The factor 2 in the above equation comes from the two-to- 
one mapping between the spin and vertex configurations 
(two spin configurations always correspond to one vertex 
configuration). 

The Boltzmann's weights, which correspond to eight 
possible line coverings of the eight-vertex model shown 
in Fig. [21 can directly be obtained from Eq. @ by sub- 
stituting a respective spin configuration of the Ising spins 

cji = [x)2 = 2 exp [ Kz H '^—j — ^ ) cosh {Kx — Ky) 



-h2exp y-K, 
UJ3 — UJ4 ~ 2 exp ( K 
2 exp ( -if. 



4 

K1+K2 
4 

K1 + K2 
4 

K1+K2 



cosh {Kx + Ky), (6) 
cosh {Kx — Ky) 



cosh (Kx + Ky), (7) 



UJ5 = uje = 2 exp ( Kz - —^—^—^^ cosh {Kx - Ky) 
+ 2exp ( -Kz + 7 } cosh{Kx + Ky), (8) 



cl'7 = ojg = 2 exp ( K 
+ 2 exp (-K, 



4 

K1-K2 



4 

K1-K2 



cosh {Kx — Ky) 



cosh{Kx + Ky). (9) 



It is quite obvious from the set of Eqs. ([1])-® that the 
Boltzmann's weights are pairwise equal to each other and 
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there are merely four independent Boltzmann's weights. 
The spin- 1/2 Ising-Heisenberg model defined through the 
Hamiltonian ^ has been accordingly mapped to the 
eight-vertex model, which quite generally satisfies the 
Baxter's zero-field condition lui — lu2, lu^ — lu^, lu^ — i-^ej 
and Loi — Us P, lUUll- Hence, it follows that the exact 
solution of the spin-1/2 Ising-Heisenberg model with the 
pair XYZ Heisenberg interaction and the quartic Ising 
interactions can be extracted from the Baxter's exact so- 
lution of the corresponding zero-field eight- vertex model 
[l], [13, HI]. For instance, the critical condition of the 
zero-field eight-vertex model 

UJi + LU3 + UJ5 + ujj = 2max{wi, W3, uj^, wy}, (10) 

directly determines phase transitions of the spin-1/2 
Ising-Heisenberg model if the effective Boltzmann's 
weights ([H])-® SLie substituted into this critical condi- 
tion. It is also worthy to mention that the critical ex- 
ponents, which characterize the phase transitions of the 
zero-field eight-vertex model, continuously change with 
the parameter /i = 2 arctan(a;5W7/a;iW3)^/^ by following 
the formulas 

, TT TT ,77 

a = a = 2 , p = , I' = I' = — , 

fi 16/i 2/z 



model and the coupling parameters i?i , R2, and i?4 of the 
zero-field eight-vertex model in the Ising representation 



7 = 7 



8^' ^ 



15, 



1 

^=4- 



(11) 



Note that the set of relations (jlip will also govern changes 
of the critical exponents of the Ising-Heisenberg model 
provided that the effective Boltzmann's weights (H])-® 
are used for a calculation of the parameter fj,. 

To provide an alternative proof of the exact mapping 
equivalence between the zero-field eight- vertex model and 
the spin-1/2 Ising-Heisenberg model, one can utilize the 
fact that the eight-vertex model on a square lattice can 
also be reformulated as two spin-1/2 Ising square lat- 
tices coupled together by means of the quartic interaction 
[24,^]. In accordance with this statement, the effective 
Boltzmann's factor ^ could be eventually replaced via 
appropriate algebraic transformation of the form 

UJp (cTpi, tTp2, CTpg, (Tp4) = 

2 exp {K^ + -ftTiCT^icr^g + if2fT^2^p4) cosh {K^ - Ky) + 
2 exp {-K^ - i^iCTpiCTps - if2crp20-p4) cosh {K^ + Ky) 

= Ro exp(i?iCrpiCrp3 + i?20-p2^p4 + ^4CTpiO-p2CTp3'^p4Xl2) 

where the mapping parameters Ri and R2 denote the 
effective pair interactions in two different Ising square 
lattices and the mapping parameter i?4 determines the 
effective quartic interaction that couples together both 
Ising square lattices. The algebraic transformation (fT^ 
must satisfy the 'self-consistency' condition, which means 
that it must hold independently of spin states of four 
Ising spins involved therein. The 'self-consistency' condi- 
tion thus provides a simple connection between the effec- 
tive Boltzmann's weights ((B])-® of the Ising-Heisenberg 



Ro^{u;ii03UJ5L07)^^^ , 
( WiWyX 



i?i=ln 
i?2 = ln 



i?4 = 41n 



(13) 
(14) 

(15) 

(16) 



This is actually an alternative proof of the exact mapping 
equivalence between the spin-1/2 Ising-Heisenberg model 
and the zero-field eight- vertex model on a square lattice. 

The usefulness of the latter equivalence consists in that 
another three useful observations can be made from it. 
First, it is quite evident from Eqs. ((M]) and (fTS]) that the 
effective pair interactions of both Ising square lattices 
become equal to each other (R\ = R2) by imposing the 
condition W5 = wy, which is equivalent to Ki — K2 (or 
Ji — J2). This means that a difference in two diagonal 
quartic Ising interactions merely causes a difference in 
the effective pair interactions of the Ising square lattices. 
Second, the quartic interaction i?4 that couples together 
two Ising square lattices vanishes only if at least one from 
both quartic Ising interactions Ki or K2 ( Ji or J2) equals 
zero provided that all interactions are finite. Another 
particular cases with the zero effective quartic interac- 
tion i?4 are the Ising and XY limit of the Heisenberg pair 
interaction. One actually obtains i?4 0, i?i — > ii^i, 
and i?2 — *■ ±if2 in the Ising limit ±cxd, whereas 

i?4 ^ 0, i?i ^ —Ki, and R2 —K2 is obtained in the 
XY limit if a: J Ky — > ±00. Apart from these rather trivial 
cases from the Ising universality class, one may expect 
that the spin-1/2 Ising-Heisenberg model will generally 
exhibit continuously varying critical exponents satisfy- 
ing the weak universality hypothesis [5§| due to the non- 
zero effective quartic interaction i?4. Third, the spin-1/2 
Ising-Heisenberg model defined through the Hamiltonian 
([1]) can in turn be refined by the terms depending on the 
Ising spin pairs (JpiOp-^ and cfp2'^pi: without disturbing the 
exact mapping equivalence to the zero-field eight-vertex 
model. As a matter of fact, an inclusion of two diagonal 
pair interactions —Ji<ypiO'p-i, —J2'^p2'^pi^ the quar- 
tic interaction — J4CpiO'p2'''p3^p4 ii^^o the Hamiltonian ([T]) 
merely adds the respective coupling constant as an aux- 
iliary constant factor into a definition of the mapping 
parameter and (fTB|) . respectively. It is note- 

worthy, however, that the model system is mapped to 
the more general and yet unsolved eight-vertex model 
by introducing the nearest-neighbor pair interactions 

-./('(cr^icr^2+^30-p4) and -J2'('7p20-p3 + f^pic^p4) into the 
Hamiltonian Even though there does not exist the 
general exact solution of the corresponding eight-vertex 
model, it is now well established that the critical expo- 
nents still vary continuously with the interaction param- 
eters even for this more complex but surely more realistic 

model [13, in, m m, m m, mm,m- 
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III. RESULTS AND DISCUSSION 

Before proceeding to a discussion of the most interest- 
ing results, it is worthy to notice that our further analy- 
sis will be restricted just to a particular example of the 
spin- 1/2 Ising-Heisenberg model with the identical quar- 
tic Ising interactions Ji — J2 ~ J4 and the more sym- 
metric XXZ Heisenberg interaction with = Jy — J A, 
Jz — J in order to avoid over-parametrization of the 
model under investigation. Furthermore, it also conve- 
nient to normalize all interaction parameters with re- 
spect to the z-component of the XXZ Heisenberg inter- 
action, which will henceforth serve as the energy unit. 
Accordingly, the dimensionless temperature will be set 
to k^T/ J, a relative strength of the quartic Ising inter- 
actions will be proportional to the ratio J4/ J, and finally, 
the parameter A = J.^/ J = Jy/ J will measure a relative 
strength of the exchange anisotropy in the XXZ Heisen- 
berg interaction. 

Let us perform first a comprehensive analysis of the 
ground state. It is worthwhile to remark that the ground- 
state spin arrangement will be thoroughly determined by 
the lowest-energy eigenstate that enters into the great- 
est Boltzmann's weight, since each Boltzmann's weight 
involves four eigenenergies that correspond to possible 
eigenstates of the Heisenberg spin pair at a given config- 
uration of the Ising spins (remember that the Ising spin 
configurations are unambiguously assigned to the Boltz- 
mann's weights according to a scheme depicted in Fig. [2]). 
It might be therefore quite useful to quote initially ex- 
plicit expressions for the effective Boltzmann's weights of 
simplified version of the spin-1/2 Ising-Heisenberg model 



oji —2 exp 
X cosh 
= 2 exp 
X cosh 



(3 (2 J + J4) 
13 (2 J - J4) 



2 exp 



2 exp 



(3 (2 J + J4) 
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(17) 



ujq — LOj — 2 exp I — ^ '^^'^ \ — 4~ ) '"'^^^ I — 2 — ) ■ 

It directly follows from the set of Eqs. ([TT]) that the 
greatest Boltzmann's weight is either loi or 0^3. Be- 
sides, another useful observation is that both these Boltz- 
mann's weights are simply connected through the rela- 
tion a;i(±J4) = a;3(=pJ4) and one may further assume 
the positive quartic Ising interaction J4 > without 
loss of generality. This simplification is a direct con- 
sequence of the affirmation that the investigated model 
system is invariant under the transformation J4 — > — J4 

and (crpi,crp2,CTp3''^p4) ^ (o'pi^ o-p2> -f^pS' -'^^4)- The 
ground-state phase diagram must be therefore symmetric 
with respect to the J4 = axis and the respective spin ar- 
rangements to emerge in a sector J4 < can be obtained 



from its symmetry related counterpart in a sector J4 > 
by a mere interchange of the Ising spin configurations 
inherent to the Boltzmann's weights loi and lo^. 

For better orientation, we will distinguish our subse- 
quent analysis of the spin-1/2 Ising-Heisenberg model 
with the ferromagnetic (J > 0) and antiferromagnetic 
(J < 0) Heisenberg pair interaction, respectively. In the 
former case, the greatest Boltzmann's weight is loi for 
A < 1 and ^3 for A > 1, under the assumptions T = 
and J4 > (note that the reverse conditions hold for the 
case J4 < 0). The global ground-state phase diagram 
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FIG. 3: Ground-state phase diagram of the spin-1/2 Ising- 
Heisenberg model with the ferromagnetic interaction J > 
in the J4/J — A plane. Each sector of the phase diagram con- 
tains typical spin configurations of the Ising and Heisenberg 
spins, which form an elementary octahedron. The sign '±' 
denotes spin states cr^ = ±1/2 and = ±1/2 of the Ising 
and Heisenberg spins, respectively. Broken lines, which con- 
nect both apical Heisenberg spins in the phases |II) and |IV), 
label the entangled spin state (|±, — ) ± | — , ±)) /\/2. 

of the spin-1/2 Ising-Heisenberg model with the ferro- 
magnetic Heisenberg interaction is displayed in Fig. [S] 
Clearly, the ground-state diagram constitute four differ- 
ent phases 



(18) 



-+))5,>(19) 
(20) 



ii) =ni+,±,+,± 

iii) =ni+,±,-,T).,-^(i+, 

p ^ 

iiii)=ni+'±'-'T)-j+,+)s, 

p 

iiv)=ni+'±'+'±)'^.75 (1+'^) + (21) 
p ^ 

which are written in a form of the product over the re- 
spective lowest-energy eigenstate of the Hamiltonian ([T]) 
that is of course identical for each elementary unit cell 
(octahedron). Note furthermore that the first (second) 
ket vector behind the product symbol determines spin 
states of the Ising (Heisenberg) spins and another equiv- 
alent representations of the eigenstates can be obtained 
from the eigenvectors (IT51) - (PT|) under the reversal of all 
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Ising and/or Heisenberg spins. The phases |I) and |III) 
are consequently eight-fold degenerate as spin reversals 
of the Ising and Heisenberg spins can be performed inde- 
pendently of each other, whereas the phases |II) and |IV) 
are just four-fold degenerate as the spin reversal of the 
Heisenberg spins does not in fact lead to a new eigenstate. 

It is worthwhile to remark that another two obvi- 
ous features directly follow from Fig. [3] and Eqs. (fT5)) - 
((2T|) . The Heisenberg spins are ferromagnetically aligned 
with respect to each other in the phases |I) and |III) 
where A < 1, while they reside the entangled spin state 
(| + , -) + |-,-f)) /\/2 in the phases |II) and |IV) where 
A > 1. In addition, there appears a perfect ferromag- 
netic or antiferromagnetic alignment on a square lattice 
of the Ising spins in the phases |I) and |IV) in contrast to 
two feasible superantiferromagnetic orderings that occur 
in the phases |II) and |III), where a perfect ferromag- 
netic arrangement of the Ising spins in a horizontal di- 
rection is accompanied with a perfect antiferromagnetic 
arrangement in a vertical direction, or vice versa. The 
overall understanding of a nature of the aforedescribed 
spin arrangements readily follows from the energy spec- 
trum of the spin-1/2 XXZ Heisenberg dimer. Namely, 
the lowest-energy eigenstate of the spin-1/2 XXZ Heisen- 
berg dimer is the ferromagnetic state |-f,+) if one con- 
siders the easy- axis exchange anisotropics A < 1, while 
the entangled spin state (| + ,— ) + | — ,+)) /v^ becomes 
the lowest-energy eigenstate for the easy-plane exchange 
anisotropics A > 1. The relevant spin arrangement of the 
Ising spins is subsequently driven by the effort to mini- 
mize the energy gain arising from the quartic Ising inter- 
actions. Thus, there must be either zero or two spins with 
opposite orientation with respect to the others among the 
four spins (two Ising and two Heisenberg) involved in the 
positive quartic Ising interaction J4 > 0, whereas there 
must be just one unaligned spin among them whenever 
the negative quartic Ising interaction J4 < is assumed. 

Surprisingly, the ground-state analysis becomes much 
more simple for the spin-1/2 Ising-Heisenberg model with 
the antiferromagnetic (J < 0) Heisenberg interaction. 
It is quite apparent from the set of Eqs. pT]) that the 
greatest Boltzmann's weight is always being uj^ {^1) if 
one considers the positive (negative) quartic Ising inter- 
action. Owing to this fact, the overall ground-state phase 
diagram constitute just two different phases 

|V)=ni+,±,-,T).,-^(| + ,-)-|-,+))5,, (22) 
p 

|VI)=ni + ,±,+,±)..^(l + ,->-|-,+»s,. (23) 

whereas the phase |V) is the ground state for J4 > 
and the phase |VI) for J4 < 0. It should be stressed 
that both four- fold degenerate phases |V) and |VI) quite 
closely resemble the phases |II) and |IV) described pre- 
viously by the analysis of the ferromagnetic model. 
As a matter of fact, the only difference between the 
phases |V) and |II) (or |VI) and |IV)) consists in a 



quantum entanglement of the Heisenberg spin pairs, 
which is described by the antisymmetric singlet-dimer 
wave function (|+,— ) — | — ,+)) /v^ in the antiferromag- 
netic model with J < and its symmetric counterpart 
(1 + ,—) + | — ,+)) /\/2 in the ferromagnetic model with 
J > 0. 

Now, let us turn our attention to a detailed analysis 
of the finite-temperature phase diagrams. It is worth- 
while to recall that phase transition lines of our simpli- 
fied version of the spin-1/2 Ising-Heisenberg model can 
be straightforwardly obtained from the critical condi- 
tion of the corresponding zero-field eight- vertex model by 
substituting the effective Boltzmann's weig hts (HZl) into 
Eq. (fTO|) . It should be also mentioned that the greatest 
Boltzmann's weight might be either wi or ^3. In the 
former case (wi > W3), the critical condition reads 



cosh ( ^^^A ) — exp ' ^'^'^ 



sinh 



- 1 



sinh 



1 



(24) 



while in the latter case {uji < lu^) the critical condition 
becomes 



(25) 



where Pc = l/(fcB7c) and Tc is the critical temperature. 
Interestingly, both critical conditions ((24|) and (|25p can 
be joined together to yield a single critical condition 



sinh 



/3c|J4| 



|x exp(Y) +cosh(Y^) 



exp 



13c J 
2 



— cosh 



(26) 



which is valid in a whole region of the parameter space. 
The absolute value of the quartic Ising interaction ap- 
pears in Eq. ((26)) as a direct consequence of the symme- 
try relation between the Boltzmann's weights a;i(±J4) — 
LU3{^J4), which are mutually interchangeable under the 
transformation J4 — s- — J4. Accordingly, the critical tem- 
perature of the spin-1/2 Ising-Heisenberg model must be 
independent of the sign of the quartic Ising interaction 
and one may still consider positive values of the quartic 
Ising interaction J4 > without loss of the generality. 
Henceforth, we will therefore restrict ourselves just to an 
analysis of the critical lines of the phases |I), |II), and |V) 
emerging for the positive quartic Ising interaction J4 > 0. 
It should be nevertheless mentioned that the critical lines 
for the phases |III), |IV), and |VI) appearing in the case 
J4 < are formally identical with the displayed critical 
lines for the phases |I), |II), and |V), respectively. 

Let us construct first the finite-temperature phase dia- 
gram of the spin-1/2 Ising-Heisenberg model with the fer- 
romagnetic Heisenberg interaction. The greatest Boltz- 
mann's weight for the ferromagnetic model with J > is 
uji if and only if 



sgn ( J4 



exp ( - cosh r^A 



> 0, 



(27) 
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otherwise ^3 becomes the greatest Boltzmann's weight. 
Consequently, the critical condition ([M)) determines 
phase transitions just if the inequality (|27|) holds, while 
the critical condition (|25p is valid in the rest of the param- 
eter space. Phase transition lines in the form of critical 
temperature versus exchange anisotropy dependences are 
drawn in Fig.[3]for three different relative strengths of the 
quartic Ising interaction. As one can see, the critical tem- 
perature generally exhibits a remarkable dependence on 
the exchange anisotropy with two marked wings of criti- 
cal lines that merge together at the ground-state bound- 
ary A = 1 between the phases |I) and |II). It should be 
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FIG. 4: Solid lines depict the critical temperature as a 
function of the exchange anisotropy for the spin-1/2 Ising- 
Heisenberg model with the ferromagnetic Heisenberg interac- 
tion J > at three different relative strengths of the quar- 
tic Ising interaction. The left wing of critical lines deter- 
mines phase boundaries for the phases |I) (J4 > 0) or |III) 
(J4 < 0), while the right wing determines phase boundaries 
for the phases |II) (J4 > 0) or |IV) (J4 < 0). The broken line 
shows the disorder solution obtained by solving Eq. (|28|) . 



the ground state and the phase |I) has slightly higher 
energy, on account of the much larger entropy gain of 
the phase |I) obtained upon the temperature increase. It 
is worthy to mention that similar re-entrant phase tran- 
sitions have been found in a variety of frustrated Ising 
models, which have been exactly solved by establishing 
a precise mapping relationship with the corresponding 
free-fermion vertex models ilfl, lZI, 72, 73, 74, 75]. To 
bring a deeper insight into the re-entrant phenomenon, 
it is also worthwhile to inspect the disorder solution to 
be derived from the condition 

cji = 0)3 <^ exp I 1 = cosh I — ^ — 1 , (28) 

where /3d = l/(fcB7b) and is the disorder tempera- 
ture. The disorder solution entails an effective reduction 
of the dimensionality and ensures a disordered nature of 
the spin system on the particular subvariety of the pa- 
rameter space given by Eq. ^ i, [ll S [tJ, [li, [ll, iB, 
pTSj . It is quite apparent from Fig. |3] that the disorder 
(broken) line calculated from the condition has a fi- 
nite positive tangent, which can be approximated in a rel- 
atively wide temperature range by a linear dependence of 
the disorder temperature To on the exchange anisotropy 
A through the relation k^T^^jJ = (A - l)/ln4. 

For the sake of comparison, the critical temperatures 
of the spin- 1/2 Ising-Heisenberg model with both ferro- 
magnetic (J > 0) as well as antiferromagnetic (J < 0) 
Heisenberg pair interaction are depicted in Fig. [5] for 
three different relative strengths of the quartic Ising in- 
teraction. The critical lines of the ferromagnetic model 
are displayed as solid lines, while the critical lines of the 
antiferromagnetic model are drawn as broken lines. The 



mentioned that the left wing of displayed phase bound- 
aries is a solution of the critical condition , while the 
right wing is a solution of the critical condition (|25p . In 
this regard, the left wings deliminate critical lines of the 
phase |I), whereas the right wings represent critical lines 
of the phase |II). The critical temperature of the phase 
1 1) at first steadily decreases as the exchange anisotropy 
A is raised from zero, then it exhibits an interesting re- 
entrant behavior before it finally tends to zero in the 
limit A ^ 1. It is noticeable that the re-entrant phase 
transitions are observable just in a relatively narrow in- 
terval of the exchange anisotropics A G (1, A^ax), whose 
upper bound Amax is the greater, the greater a relative 
strength of the quartic Ising interaction is. Contrary to 
this, the right wing of critical lines, which allocate phase 
transitions of the phase |II), characterizes a monotonous 
increase of the critical temperature with the exchange 
anisotropy. An origin of the observed re-entrant tran- 
sitions obviously lies in two times higher degeneracy of 
the phase |I) compared to the phase |II). Accordingly, 
two successive re-entrant transitions from the paramag- 
netic phase to the phase |I) and back may take place in 
the parameter region, where the phase |II) constitutes 



0.3 




A 

FIG. 5: The critical temperature as a function of the ex- 
change anisotropy for the spin-1/2 Ising-Heisenberg model 
with the ferromagnetic (solid lines) or antiferromagnetic (bro- 
ken lines) Heisenberg pair interaction at three different rela- 
tive strengths of the quartic Ising interaction. 

displayed critical lines of the antiferromagnetic model in 
fact determine critical temperatures of the phase |V), 
which is the only possible ground state when J < and 
J4 > 0. Evidently, the antiferromagnetic model generally 
exhibits a rather small variation of the critical tempera- 
ture upon varying the exchange anisotropy in comparison 
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with the marked two- wing dependence of the critical tem- 
perature that shows the ferromagnetic model. Another 
interesting fact to observe here, and also derived from 
Eq. ([^5]) . is that the critical temperatures of the ferro- 
magnetic and antiferromagnetic model become equal to 
each other in two limiting cases A ^ and A — > oo. 
From this perspective, the most obvious difference be- 
tween the critical temperatures of the ferromagnetic and 
antiferromagnetic Ising-Heisenberg models appears in a 
vicinity of the isotropic limit A = 1, where the critical 
temperature of the ferromagnetic model asymptotically 
vanishes due to the zero-temperature phase transition be- 
tween the phases |I) and |II). 

At this point, let us make certain comments on changes 
of critical exponents along the critical lines. For this 
purpose, typical changes of the critical exponent a along 
the phase transition line of the spin- 1/2 Ising-Heisenberg 
model with the ferromagnetic Heisenberg interaction and 
I J4I/J = 1 are depicted in Fig. [S] Sohd lines display a 
variation of the critical temperature with the exchange 
anisotropy, which are scaled with respect to left axes, 
while broken lines scaled with respect to right axes show 
the relevant changes of the critical exponent a. As could 
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FIG. 6: The changes of the critical exponent a along the 
critical line of the spin-1/2 Ising-Heisenberg model with the 
ferromagnetic interaction J > and | J4I/J = 1.0. Solid lines, 
which are scaled with respect to left axes, show the critical 
temperature as a function of the exchange anisotropy. Bro- 
ken lines, which are scaled with respect to right axes, display 
the relevant changes of the critical exponent —a in a semi- 
logarithmic scale. Fig. |6Jd) illustrates a detail of the parame- 
ter region, where both re-entrant phase transitions as well as 
a quantum critical point occur. 

be expected, the critical exponent a varies continuously 
along the line of critical points according to the rela- 
tions (fTT|) . which connect changes of the critical expo- 
nents to respective changes of the interaction parameters 
through the parameter fi. It turns out that the inves- 
tigated model system generally exhibits rather smooth 
continuous phase transitions, since the order of phase 
transition is proportional to r = 2 — a (see for instance 
pp. 16-17 in Ref. [l|) and the critical exponent a lies 
within the range a € (—00, 0). The most striking finding 
is, however, that the critical exponent a exhibits a pecu- 
liar singularity in a vicinity of the ground-state boundary 



between the phases |I) and |II) where a ^ — 00. Owing 
to this fact, the zero-temperature phase transition be- 
tween the phases |I) and |II) emerging at A = 1 is of 
infinite order and hence, this special critical point actu- 
ally represents a remarkable quantum critical point. It is 
also worthy to notice that the displayed variations of the 
critical exponent —a quite closely resemble dependences 
of other critical exponents /?, 7 and v due to a similar 
mathematical structure of the relations pTjl . 

Furthermore, the changes of the critical exponent a 
along the critical line of the spin-1/2 Ising-Heisenberg 
model with the antiferromagnetic Heisenberg interaction 
are shown in Fig. [T] In this particular case, the changes 
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FIG. 7: The changes of the critical exponent a along the crit- 
ical line of the spin-1/2 Ising-Heisenberg model with the anti- 
ferromagnetic interaction J < and |J4|/|J| = 1.0. The solid 
line, which is scaled with respect to the left axis, shows the 
critical temperature as a function of the exchange anisotropy. 
The broken line, which is scaled with respect to the right axis, 
displays the relevant changes of the critical exponent a along 
this critical line. 

of the critical exponent a are restricted to a rather nar- 
row finite interval even though a negative value of the 
critical exponent a still implies a continuous nature of 
the phase transitions r > 2. It is also quite interest- 
ing to notice that the critical exponent a monotonically 
increases from its smallest value at A = before it grad- 
ually tends towards a ^ in the limit A — > 00. Al- 
together, the critical lines of both ferromagnetic as well 
as antiferromagnetic model turned out to be the lines 
of continuous phase transitions of the order r > 2, but 
the ferromagnetic model generally exhibits much more 
pronounced changes of the critical exponents including a 
rather strange quantum critical point of the infinite order 
with diverging critical exponents. 

Before concluding, let us make few remarks on a crit- 
ical behavior of another particular case of the spin-1/2 
Ising-Heisenberg model with Jx = Jy = </A, Jz = J and 
the quartic Ising interactions Ji — —J2 = J^- This par- 
ticular case differs from the previous one just in a nature 
of the quartic Ising interactions, which are of equal rela- 
tive strengths but of different signs. In this respect, the 
positive quartic interaction will prefer spin alignments 
with either zero or two opposite spins among the four 
spins (two Ising and two Heisenberg) involved therein, 
while the negative quartic interaction will favour spin 
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alignments with just one unaligned spin among them. 
Interestingly, the phase diagrams of this particular case 
will be essentially identical with the ones discussed pre- 
viously, since there is a simple relation between the effec- 
tive Boltzmann's weights of both particular cases. As 
a matter of fact, it directly follows from Eqs. ([SI)-® 
that the effective Boltzmann's weights of the model with 
two non- uniform quartic interactions Ji = — J2 = J4 
can be expressed through the Boltzmann's weights (fT7|) 
of the uniform case as uj[ = ^3 = oj^, uji^ = W3, and 
Lu''^ = uji. It is quite apparent that the role of differ- 
ent Boltzmann's weights is merely interchanged and the 
greatest Boltzmann's weight is now being either or 
Wg. However, this fact does not affect the ground-state 
and finite-temperature phase diagrams because the crit- 
ical condition (jlOp is quite symmetric with respect to 
all four Boltzmann's weights involved therein. Hence, it 
follows that the zero- and finite-temperature phase dia- 
grams shown in Figs. [3][5] remain valid and one should 
only perform a respective change of the Ising spin config- 
urations from wi — > uj'.j and W3 —fco'^, respectively. The 
ground-state phase diagram of the ferromagnetic (J > 0) 
model with the non-uniform quartic interactions, which 
is shown in Fig. [3l thus constitute the phases 

|l'>=ni + 'T,+,±>.J+,+>5„ (29) 
p 

|Il')=n l±' +' -> + (30) 

p ^ 

iiii'>=ni±.+.T,+>.j+,+>5„ (31) 

p 

|IV') =n 1+' T, +, ±).,^ (I+, -) + I-, +))s,- (32) 
p ^ 

while the ground state of the antiferromagnetic (J < 0) 
model is either being |V') for J4 > or |VI') for J4 < 

IV') =n T, ^ (|+, -) - I-, +))s,, (33) 
p ^ 

|VI') =n I+' +, ±)., ^ -) - +»5,- (34) 

Despite similarity of the ground-state and finite- 
temperature phase diagrams, the critical exponents of 
the particular case with the non-uniform quartic inter- 
actions are fundamentally different from the ones of the 
former particular case with two identical quartic inter- 
actions. Namely, the changes of the critical exponents 
depend through the set of Eqs. pT|) on the parameter /i, 
which is not symmetric with respect to all four Boltz- 
mann's weights. Therefore, it might be quite interest- 
ing to ascertain how the critical exponents change by 
assuming two quartic interactions of different sign (na- 
ture). The variations of the critical exponent a along 
the critical line of the ferromagnetic model (J > 0) are 
displayed in Fig. [5] for the particular case with | J4I/ J = 1. 
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FIG. 8: The changes of tlie critical exponent a along the 
critical line of the spin- 1/2 Ising-Heisenberg model with the 
ferromagnetic interaction J > and the relative strength of 
the non-uniform quartic interactions j J4I/J = 1.0. Solid lines, 
which are scaled with respect to left axes, show the critical 
temperature as a function of the exchange anisotropy. Broken 
line, which are scaled with respect to right axes, display the 
relevant changes of the critical exponent a. 



It is quite obvious from this figure that the critical ex- 
ponent a exhibits an outstanding dependence with the 
global maximum ckmax = 1 emerging at A = 1. As a re- 
sult, the zero-temperature phase transition between the 
phases |I') and |II') (or |III') and |IV')) is in fact a dis- 
continuous first-order phase transition [r — 1) unlike the 
aforedescribed continuous phase transition of the infinite 
order between the phases |I) and |II) (or |III) and |IV)). 
Besides, the critical exponent a evidently acquires posi- 
tive values from the interval a G (0, 1), which indicates 
discontinuous nature of the phase transitions (r < 2) 
along the whole critical line of the ferromagnetic model 
with the non-uniform quartic interactions. 

For completeness, we depict in Fig. [9] the critical ex- 
ponent a of the antiferromagnetic model (J < 0) with 
the same relative strength of non-uniform quartic inter- 
actions |J^4|/| J| = 1. It can be clearly seen from this fig- 
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FIG. 9: The changes of the critical exponent a along the crit- 
ical line of the spin- 1/2 Ising-Heisenberg model with the an- 
tiferromagnetic interaction J < and the relative strength of 
the non- uniform quartic interactions | J4I/I J| = 1.0. The solid 
line, which is scaled with respect to the left axis, shows the 
critical temperature as a function of the exchange anisotropy. 
The broken line, which is scaled with respect to the right axis, 
displays the relevant changes of the critical exponent a. 

ure that the critical exponent a monotonically decreases 
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with the increase of the exchange anisotropy even al- 
though it always remains positive. It is noteworthy that 
the precisely opposite trends were observed in the for- 
mer particular case with the uniform quartic interactions 
and there are only two common features of both par- 
ticular cases with the antiferromagnetic pair interaction 
J < 0. First, the changes of the critical exponent a are 
restricted just to a rather narrow finite interval and sec- 
ond, the same asymptotic value a = is achieved in the 
limit A ^ oo. Altogether, it could be concluded that 
the Ising-Heisenberg model with the non-uniform quar- 
tic interactions generally exhibits a discontinuous phase 
transitions of the order r < 2 no matter whether the fer- 
romagnetic or antiferromagnetic nature of the Heisenberg 
pair interaction J is assumed. 

IV. CONCLUSIONS 

In the present work, we have furnished proof of an 
exact mapping equivalence between the spin-1/2 Ising- 
Heisenberg model on a two-dimensional lattice of edge- 
sharing octahedrons and the zero-field eight-vertex model 
on a square lattice. In accordance with this exact map- 
ping correspondence, the critical behavior of the model 
under investigation closely resembles an outstanding crit- 
ical behavior of the zero-field eight-vertex model and 
Ashkin- Teller model [l|, which exhibit critical lines with 
continuously varying critical exponents satisfying the 
weak universality hypothesis [59|. It is worthwhile to 
remark, moreover, that the similar precise mapping rela- 
tionship between the spin-1/2 Ising-Heisenberg model on 
a square-hexagon lattice and the zero-field eight-vertex 
model was recently found by Valverde et al. in a re- 
stricted region of the interaction parameters [t^ . Unlike 
this case, the exact mapping correspondence with the 
zero-field eight-vertex model holds in our case quite gen- 
erally, i.e. it is not restricted to any particular subvariety 
of the parameter space. 

Exact results for the spin-1/2 Ising-Heisenberg model 
with the pair XYZ Heisenberg and quartic Ising inter- 
actions imply that the model with the antiferromag- 
netic pair interaction surprisingly exhibits less significant 
changes of both critical temperatures as well as critical 
exponents than the model with the ferromagnetic pair 
interaction. The most interesting finding to emerge from 
the present study is however an exact evidence of a quan- 
tum critical point of the infinite order, which character- 
izes a peculiar singular behavior of the critical exponents 
in a close vicinity of the isotropic limit of the Heisenberg 
pair interaction. In addition, it was shown that the crit- 
ical exponents vary continuously over the entire range of 
allowed values by changing the exchange anisotropy in 



the Heisenberg pair interaction and the relative strength 
of the quartic interaction. From this point of view, the 
investigated Ising-Heisenberg model represents a rare ex- 
ample of the exactly solved quantum spin model with 
such an unusual weak-universal critical behavior. 

Next, it is worth noticing that the Ising-Heisenberg 
model with the pair Heisenberg and quadratic Ising in- 
teractions can also be interpreted as an interesting ex- 
ample of the exactly solved quantum dimer system with 
a rather complicated even-body interactions between 
dimers. This means that the system of quantum dimers 
in some specific fiuctuations or fields is also equivalent 
with the exactly solved zero-field eight-vertex model, 
which contradicts the standard universality hypothesis in 
that its critical exponents vary continuously with the in- 
teraction parameters. In addition, the quartic and other 
even-body interactions turned out to play an important 
role in determining magnetic properties of several insu- 
lating magnetic materials [3ij 32, 33, 34, i55, 36, 37, 3^, 
which makes the presented exact results more interesting 
also from the experimental point of view. Even though 
it would be rather striking coincidence if some real mag- 
netic material would obey very specific topological re- 
quirements of the model under investigation, it is quite 
reasonable to suspect that our exact results might at least 
shed light on some important aspects of the critical be- 
havior of real magnetic materials. 

Last but not least, it should be mentioned that it 
could be quite interesting to explore also temperature 
variations of some basic thermodynamic quantities (such 
as entropy, specific heat, etc.) that might exhibit a 
rather spectacular thermal dependences especially near 
the zero-temperature phase transitions. Furthermore, 
several interesting extensions and generalizations of the 
present version of the Ising-Heisenberg model come into 
question besides the most obvious ones mentioned at the 
end of Section [Hi For instance, it is possible to extend 
the present Ising-Heisenberg model by including higher- 
order triplet, quintuplet, and sextuplet interactions be- 
tween the Ising and Heisenberg spins, or to solve exactly 
the analogous Ising-Heisenberg model with the Heisen- 
berg spins S > 1/2 that accounts for another interaction 
terms like the single-ion anisotropy and the biquadratic 
interaction. In this direction will continue our next work. 
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